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Abstract 

We exploit the local loop dynamics calculated in prepotential formulation to compute the 
pertrubation expansion in the strong coupling limit of lattice gauge theory. A new exact simulation 
technique is developed to simulate all possible loop states on an infinite lattice using exact BFACF 
algorithm. This loop perturbation numerical calculations are free from any limitation on lattice 
size. Using this, we calculate the correction of string tension for SU(2) theory in 2+1 dimension 
upto order g 64 . 


1 Introduction 

Lattice gauge theory, from its beginning has been extremely popular and advantageous primarily 
because of the plausibility of practical computations of high energy physics which was almost impossible 
to address in continuum. Just after the advent of lattice gauge theory in its Euclidean version pQ, 
a Hamiltonian version was also proposed in [2]. However in the course of next three and a half 
decades, the later remains much unexplored primarily because of its not so well adaptation in numerical 
simulations, at least upto the level in which the first one has achieved. But, both the Hamiltonian and 
Euclidean path integral approach of lattice gauge theory has their own benefits over the other. The 
Hamiltonian formulation allows one to directly work with the fock space, vacuum structure and the 
spectrum of the theory which is not such straight forward in the other way. 

Lattice gauge theories do suffer from a huge redundant gauge degrees of freedom in general which 
is a primary obstacle in terms of computation cost. It is desirable to work within the physical degrees 
of freedom of the theory which is gauge independent or invariant. One can always achieve that either 
by making a suitable gauge choice or by formulating the theory completely in terms of gauge invariant 
variables such as Wilson loops and strings. The later prescription, although attempted by many 
theorists in the past [3], has not been of much practical advantage as the loops formulation itself 
suffers from severe non-locality and proliferation of loop and string space. 
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There has been a recent approach in the Hamiltonian lattice gauge theory in terms of local loop 
degrees of freedom 0016]. This formulation is named as prepotential formulation, where prepotentials 
are basically harmonic oscillators belonging to the fundamental representation of the gauge group and 
are defined at each lattice site. These oscillators act as the fundamental operators of the theory in terms 
of which the original canonical variables are reconstructed. By definition, prepotentials are located at 
each lattice site. The most important feature of this formulation is that, in terms of prepotentials, the 
non-abelian part of the gauge group becomes ultra local, confined to each lattice site resulting local 
loop space or gauge invariant Hilbert space of the theory situated at each site. However, the original 
non-local loops such as Wilson loops and strings do exist in the new formulation and are constructed 
by abelian weaving of local loops along neighbouring lattice sites. Prepotential formulation enables one 
to construct exact orthonormal loop basis locally at each site of the lattice and compute the dynamics 
of the theory site by site. 

It is worth mentioning that, there has been only a limited progress towards the numerical investiga¬ 
tions in Hamiltonian lattice gauge theory since the advent of the subject compared to that in Euclidean 
version of the theory. Most of the numerical studies within Hamiltonian framework till date uses the 
variational approach |7]. However, the variational study primarily depend on the trial wave function 
which is not chosen systematically such that a systematic improvement in the understanding of the 
theory can be obtained. Some progresses have been done by applying many body techniques such as 
cluster expansion, specifically coupled cluster or linked cluster algorithms [8] for lattice gauge theories. 
These cluster algorithms however involve abrupt truncation of the Hilbert space of the theory. Yet, 
they are advantageous in the sense that, they are in terms of gauge invariant variables or loops and 
hence free from any gauge redundancy. Note that, although classical Monte Carlo simulations provide 
a very powerful and accurate method for the study of Euclidean lattice gauge theories, the quantum 
Monte Carlo simulations |9] have been explored to a limited extend in the context of Hamiltonian 
lattice gauge theory. Moreover, these studies involve gauge field configurations likewise the classical 
Monte-Carlo for euclidean version of the theory and it is not possible to enforce Gauss’s law explic¬ 
itly in these calculations. Again, Monte-Carlo simulations in Hamiltonian framework, is not suitable 
to compute the string tensions and mass gaps directly as Hamiltonian eigenvalues, as one does in the 
strong-coupling representation. But in this approach, one is forced back to the more laborious approach 
used in Euclidean calculations by measuring the appropriate correlation functions and estimating the 
mass gap as the inverse of the correlation length. 

In this scenario, it is expected that, a better choice of strong-coupling basis, such as the ‘loop 
representation should avoid these problems and provide with a better numerical simulation scheme to 
explore Hamiltonian lattice gauge theory, which has not yet been demonstrated in literature so far. In 
this work, we proceed in this direction and develop a numerical simulation in terms of local orthogonal 
loop Hilbert space, whose dynamics is given in prepotential framework. 

This work is the very first attempt to exploit the novel features of prepotential formulation in 
terms of explicit computations. This local loop formulation of lattice gauge theory do exist for any 
arbitrary SU(N) gauge group. But in this work we will concentrate on the simplest non-trivial case of 
2 + 1 dimensional non-abelian pure SU(2) gauge theory defined on lattice. We will take a semi-analytic 
approach in doing so. We simulate all possible loop configurations of the theory, starting from the 
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strong coupling ground state numerically using a new exact simulation technique developed by us 
and calculate the Rayleigh-Schrodinger perturbation expansion upto a high enough order (correcting 
upto g ~ 64 ) using the analytically calculated dynamics of the theory in prepotential formulation. This 
scheme of calculation is completely gauge invariant and is defined on a truly infinite spatial lattice. 
This saves the computation cost to a huge extent. Note that, the work presented in this paper is 
simulated and calculated on an ordinary PC within a few minutes. 

The organization of the paper is as follows: We start with a brief outline of the loop perturbation 
expansion in Hamiltonian lattice gauge theory, which is Rayleigh Schrodinger perturbation expansion 
in section 2. In section 3, we introduce prepotentials and illustrate the Hamiltonian dynamics in 
terms of prepotential which is local. In section 4, we develop the algorithm for simulation of all 
loop configuration and illustrate how to construct the network of states, which is to be exploited in 
calculation of strong coupling perturbation expansion for the theory. We also illustrate the novelty 
of this algorithm by computing axial string tension semi-analytically using our algorithm. Next, in 
section 5, we discuss the results and finally conclude. 


2 Loop Perturbation Expansion in Hamiltonian Lattice Gauge 
Theory 


In this work we approach the Hamiltonian lattice gauge theory in its strong coupling limit in a complete 
gauge invariant way by numerical simulations of gauge invariant configurations of the theory. For the 
sake of completeness, we briefly discuss the Hamiltonian for the theory and the perturbation corrections 
to its spectrum in strong coupling limit. Note that, in this work we will confine ourselves to the SU(2) 
gauge group and 2 + 1 dimension for simplicity. In 2+1 dimensions, we consider discretization in the 
two spatial dimension only, keeping time continuous. The lattice is just a 2-dimensional square lattice. 

The Lattice gauge theory Hamiltonian [2 is basically a many body Hamiltonian constructed in 
terms of the canonical conjugate variables, i.e the left or right color electric field E^i R (n,i) and link 
operator C/(n,i) located on each and every link (n,i) originating at site n along direction i of the 
lattice, given by 

H — 9 Elinks — ^2 (^Tr Uplaquette + Tr ^plaquette) W 

links plaquettes 

In (|TJ) , kpi a q Ue tte = i)U(n + j)W (n + j,i)W(n : j) is product over links around a plaquette. 
The canonical conjugate variables of the theory satisfy the commutation relation: 


In ([ 2 ]), a(= 1,2,3) is the color index for SU(2) and are the generators in the fundamental rep¬ 
resentation of SU( 2) and satisfy: [^-, ^-] = ie abc ^-. The left and right electric fields follow SU(2) 


[El(n,i),U a fi (n,i)] 


U (n, i) 


[E^(n + i,i),Up(n,i)\ = 


U(n,i) y 


( 2 ) 
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algebra: 

[El(n,i),E h L (n,i)\ = ie ahc E c L {n,i), 

*)> *)] = i^abcEji{n,i), (3) 

[El(n,i),E h R (m,j)] = 0. 

Note that the right generators E^{n + i,i) are the parallel transport of the left generator E^(n,i) on 
the link (n,i) as Er(ji + i,i) = —U^(n,i)EL(n,i)U(n,i), implying, 

3 3 3 

dunks e E W", 0 = E *£(">*)*£(«, 0 = E E n( n + ^ *)£«(« + *. 0- (4) 

a=l a=l a=l 

which appear in 0 . 

In the strong coupling limit of the theory, where g 2 —» oo, the Hamiltonian is dominated by the 
electric part, i.e 

He=g 2 Y. E Unks ( 5 ) 

links 

The prime advantage of formulating gauge theories on lattice is that it is naturally confining at the 
strong coupling limit. One can perform the strong coupling perturbation expansion around this limit. 
In Hamiltonian formulation it is natural to carry out Rayleigh Schrodinger perturbation series expan¬ 
sion for lattice gauge theories with the unperturbed Hamiltonian 

H 0 =H e =g* E dunks (6) 

links 

and interaction Hamiltonian given by 

Hi = Hmag = — ^ (TXUpiaquette + ^^llaquette) ’ (V 

plaquettes 

At this point, we scale the Hamiltonian in the following way: 

H^g 2 H=J2 E Lks ~j, E (™ plaquette + 1 W} laquette ) mH 0 - (8) 

links plaquettes 

V ---' - v -- 

Ho Hl 

The unperturbed Hamiltonian can be solved exactly, with the gauge invariant electric flux states or 
loops (or string) states as eigenstates and eigenvalues counting the net flux flowing around that loop. 
This sector of the theory is completely confining showing finite mass gap and constant value of string 
tension. 

Let us consider the unperturbed Hamiltonian and its spectrum first. Let \rii) denote the i th excited 
eigenstate of Hq with energy eigenvalue Eq, where the 0 index refers to unperturbed energy, i.e 

Holm) =EZ i \n i ), for i = 0,1, 2, 3,. (9) 

for i = 0, we denote the strong coupling vacuum states by |0). The complete loop space is built 
by the successive action of the interaction Hamiltonian on |0). The perturbation correction to the 
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unperturbed spectrum is obtained from the full spectrum. For strong coupling lattice gauge theory, 
all odd order corrections do vanish and the net correction is obtained as 



( 10 ) 


The even (2 M) order corrections are obtained by the expression: 



(ii) 


In the next section we discuss and evaluate this perturbation expansion in prepotential formulation, 
where we can calculate the dynamics of any arbitrary loop state locally at each site. 


3 Prepotentials and Hamiltonian Dynamics 


The recently developed prepotential formulation of Hamiltonian lattice gauge theory has 


introduced a new level of simplicity in the theory both in terms of theoretical understanding and 
practical computations. In this novel formulation, all the gauge invariant states or loop states has a 
local description forming an orthogonal and complete basis. Moreover, the dynamics is also local. 

In this work, we will exploit the prepotential formulation to calculate the strong coupling pertur¬ 
bation expansion for pure lattice gauge theory. We approach the problem semi-analytically, or to be 
more specific, our approach consists of two steps. First, we compute the local dynamics of loop states 
analytically and then numerically simulate all possible gauge invariant configurations and perform the 
sum of their contribution to the perturbation expansion numerically. The reason for going beyond 
analytic treatment is the fact that the number of loop configurations diverge exponentially (which is 
clearly shown in our simulation) as one attempts to make higher and higher order corrections. Cal¬ 
culation of the degeneracy factor at each excited level is itself an analytically-unsolved problem in 
combinatorics. However we propose an algorithm, where all loop configurations are generated starting 
from strong coupling vacuum. In this particular work, for the purpose of illustration of the usefulness 
of the algorithm to compute physical quantities, we consider the axial string tension to be calculated 
from the string state connecting a static quark and anti-quark pair as the physical loop configuration. 

Before entering into the dynamics of the specific system, it is worth having a brief review of the 
prepotential formulation for SU(2) lattice gauge theory in 2 + 1 dimension, which sets the platform of 
this particular work. 

SU(2) prepotentials are basically two sets of Harmonic oscillator doublets (L/i?), a = 1,2 
associated to both the ends ( L/R ) of all the links on the lattice. The canonical conjugate variables in 
the Kogut -Susskind Hamiltonian, namely, the color electric fields and link operators on a link (n, i) 
are rewritten in terms of the prepotentials in the following way, keeping all the commutation relations 
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invariant: 


Left electric fields: 

El 3 +L)ya(L), 

(12) 

Right electric fields: 

tn 

III 

§ 


Link operator: 

U% = -^=fa^(L)al(R) + a a (L)a 0 (R))^= 

\Jn +1 V M / \Jn + 1 

(13) 


The above link operators are defined to satisfy the unitarity condition given by 

UU ] = U ] U = 1 & Det U = 1 (14) 

This novel prepotential approach has the following properties: 

• Each link is subdivided into a left part and a right part. Even the link operator, in matrix form 
can also be written as 0 u = U L U R . 

• The only connection between the left and right part of a link is via the abelian constraint 

n(L) = a f (L) ■ a(L) = h(R) = a f (R) • a(R ) = h (15) 


derived from the condition E\ = E^, 


In summary, the prepotential formulation makes the non abelian part of the gauge group confined 
locally at each site and abelian fluxes flowing across neighbouring sites. As a result, this approach 
contains an ultra local description of loop states, operators and their dynamics as discussed in detail 
in 000. In context of this particular work, we highlight a few important observations regarding the 
dynamics of the theory as we are interested in the case of SU(2) gauge group and 2 + 1 dimension. 


It is clear from the redefinition of link operator in terms of prepotentials in (13), that each link 
operator is a sum of two terms. One of the terms contains only prepotential creation operators causing 
net increase the flux along the link by one unit and the other contains only prepotential annihilation 
operators which decrease the flux along that link by one unit, pictorially described in figure [l] The 


+ 


Figure 1: Pictorial representation of link operator in prepotential formulation 

magnetic part of the Hamiltonian is defined around each plaquette consisting of four such link operators. 
Hence the full plaquette operator for a square lattice consists of 2 4 = 16 terms as shown in figure [ 2 ] 

We now discuss the dynamics of loop states under this plaquette operator. The loop states in 
prepotential formulation are characterized by linking numbers at each site [6|. At each site of a 
two dimensional spatial lattice, there exist six distinct linking numbers in order to characterize any 
arbitrary loop state together with three different constraints yielding three physical loop degree of 
freedom per site, as expected. The detailed dynamics of such orthonormal loops has been discussed in 
great detail in [6]. The dynamics around a plaquette can be realized more geometrically in terms of 
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Figure 2: The magnetic Hamiltonian for SU(2) lattice gauge theory 


fusion variables as described in [6] which is again shown to be exactly equivalent to the same obtained 
in terms of original link variables. As we confine ourselves to gauge group SU(2) in this work, we do 
not enter into any of these complications but just quote the required dynamics discussed in [4]. The 
complete and orthonormal loop basis for this case is characterized by 3 angular momentum quantum 
numbers per lattice site given by |ji, ji 2 )and hence, around a plaquette ‘abed’ as 


I jabcd) — h'l ? 3% j 3\ > 4*2 > 3l2) X 1 3l > 32 5 3 ±5 4*2 ’ ^' 12 ) X \jl > J2i 3\ > 3%’ ^' 12 ) X \jl j 32 j 3\ ’ $2 ’ 3l2) 


•d -d -d -d -d 


We can further identify 


3i = j\ = 3 1 
3l = 3i = 31 


32 = if = 32 
il = 32 = 3 2 


(16) 

(17) 


as the four fluxes around a plaquette. The dynamics of such states under the plaquette action are 
obtained as [4] 


{jabcd\T'rU a l )C( l \ jabcd) — -M-abcd 



312 Jl2 


3 1 


Ad Ad 

3 12 3 12 


Jl 



il2 312 


3 2 3i 


Ad Ad 

3 12 J12 


.7 2 


1 

2 

Ad 

d\ 


(18) 


In Alaftcd = D abcd N abcd Pabcd where: 


Dabcd djs , jf <% 2 , jf 2 J jf % Sjcjc Sjc 2 , ; c 2 SjdjdSjdjd , 

Nabcd H (jl; Jl; J 2 ; J 25 J3 5 Jl 5 J 2 ; J* 2 ; Jl 2 > «?125 «? 12 > J 12 ) 

Pabcd = H-l) il+h+ibl+i h-l) jl+h+jt+i *&GlJ^ 

A(il2> Jl2; 2 -?12) p- 


(19) 


In (19), D a i) C d describes the trivial <5 functions over the angular momenta which do not change under 


the action of the plaquette operator, Alfred and P a ^cd give the corresponding numerical and the phase 
factors respectively. The multiplicity factors are: n(x,y,...) = (2x + l)(2y + 1)... and A (x,y,z) 

represent the phase factors associated with a triangle with sides x, y, z: A (x,y,z) = (—1 y+v+ z =>* 
A(x,y,z) = ±1. 

Having set the dynamics of Hamiltonian lattice gauge theory in terms of the prepotentials, let us 
now concentrate on the particular problem that we are interested in for this work. In order to calculate 
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the string tension, our aim is to calculate the ground state energy of a string state connecting a static 
quark anti-quark pair as shown in figure [3] 







Figure 3: A string state connecting a static g, q pair 


To make a gauge invariant configuration, the flux string staring at q and ending at q must carry 
one unit of flux, denoted by | j = 1/2). The unperturbed energy eigenvalue for such a state is obtained 
in prepotential formulation as, 


^ links 

links 


1-7 V^)along the string 


- x length of the string in lattice units 


( 20 ) 


Hence, the lowest energy configuration of such a string state is the string connecting the quark- 
antiquark pair of minimum length, i.e along the shortest distance of the two points where the static 
sources are residing. The interaction Hamiltonian will act on such a string state and will deform it 
to another string state, which is again eigenstate of the unperturbed Hamiltonian with certain eigen 
value depending upon its length. 


At this point we choose an effective Hamiltonian from the original Kogut-Susskind Hamiltonian for 
which the interaction part is denoted pictorially in figure [4j Such an effective interaction Hamiltonian, 



plaquettes 




I 

I 

-I 


Figure 4: Effective interaction Hamiltonian 


when acts on an arbitrary string state, replaces the links on the string by a staple or vice-versa. Hence, 
the resultant string state and the original string state do always have a difference in their length by 
two units and belong to two consecutive excitation level. Such an action is illustrated in figure [5] and 
[6j Clearly, the state with the minimum length of string, say Zq is the ground states of the system with 


q 


t 





Figure 5: An arbitrary string state 


g #_n — 1-1 - *q\ V 

Figure 6: A new string state differing by a plaquette 
action 


unperturbed energy Eq = 3Zo/4. The first excited states are of length Zo + 2 with energy Eq = |(Zo + 2). 
In general n th excited state has its energy eigenvalue Eft = |(Zq + 2n). The chosen effective interaction 
Hamiltonian, acting on a particular string state at a particular unperturbed energy level moves it to 
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the next lower or upper energy states respectively. However, as the interaction Hamiltonian is summed 
over all plaquettes of the lattice, while acting on a string state, it actually create linear combination of 
many different states, all of which has the same unperturbed energy. Hence, even if the ground state 
is non-degenerate, the degeneracy increase with increasing energy leveQ 

Let us now concentrate on the matrix element of the interaction Hamiltonian between two string 
states at two adjacent energy levels. This expression is calculated analytically within prepotential 
framework. But as the number of string states at each excited level grows rapidly, one needs to 
incorporate numerical simulations of all such configurations finally yielding a semi analytic treatment 
of the problem. 


As stated earlier, the dynamics in prepotential formulation is calculated locally at each site as 


given in the expression (18). In the perturbation expansion, we have the terms (rii\Hi\nj). The states 
| rii) and | rij) are identical except around one single plaquette touching the string (an example of such 
states are shown in figure [5] and [6| . Whenever these two states are identical, they will contribute a 
factor of 1 in the matrix element. The non-trivial contribution, (around only one plaquette around 
which these two states actually differ at each step) calculated using ( [l8| ) yields a factor of 1/2. The 
dynamics obtained in terms of 1 8j symbol is rotationally symmetric and hence yield the same factor 
for one plaquette deformation of the string along any direction. 


Hence we have, for any two string states | rii) and | rij), 


{ n i\H I \nj) 


- if the two states differ by only one plaquette 
0 otherwise 


( 21 ) 

( 22 ) 


In the next section we prescribe a new numerical simulation technique to simulate all possible loop 
states of lattice gauge theory. Although, in this paper we will confine ourselves to the above discussed 
string states connecting static quark-antiquark pair, the general technique applies to simulation of all 
possible loop configurations of the theory. 


4 Numerical Simulation of Loop States 

Let us now prescribe a new simulation technique for simulating all possible gauge invariant configura¬ 
tions on the lattice starting from the strong coupling vacuum. The action of the interaction Hamilto¬ 
nian on the string (or loop) states can be numerically implemented by BFACF (Berg-Forster-Aragao- 
Caracciolo-Frohlich) algorithm [TO] . The fundamental BFACF moves were originally constructed in 
context of knotted polygons. We observe that, the entire lattice loop dynamics can indeed be imple¬ 
mented numerically by BFACF moves starting from any arbitrary loop states. However, in this work, 
we have chosen a reduced effective Hamiltonian, yielding the dynamics simulated by only the length 
changing BFACF moves for any loop configuration. Keeping the fixed length BFACF moves as well 
would have simulated the dynamics of the original Hamiltonian instead of the chosen effective one. 
Precisely, BFACF algorithm consists of a set of moves, which makes the deviation from one possible 
path connecting any two points (or even a closed path) on a lattice to another. The path can be 
1 We will show in the next section, that in reality, this degeneracy factor rises exponentially with excitation. 
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considered as a set of links connected with each other as well. Now, the length changing BFACF 
moves allow one particular link to change to a set of three links as shown in figures Mira and 
the rotations such that the length of that particular path changes by two units. In this algorithm, 
depending upon the relative orientation of three consecutive links, the middle one gets modified as 
shown in the figure. Note that, application of these moves on each and every link of a particular loop 
exactly simulates each and every loop configurations at the next excited level, and hence using this 
exact algorithm one can produce each and every configurations relevant for the problem. Being an 
exact simulation technique, it is free from any statistical error as well. 



~LT 





Figure 7: BFACF Move: 1 


Figure 8: BFACF Move: 2 



Figure 9: BFACF Move: 3 


Figure 10: BFACF Move: 4 


In this work, the configuration of interest is the string state connecting a quark anti-quark pair. 
We choose our lattice axis in such a way, that both the sources he on the X-axis. Hence, the shortest 
distance between them is always a straight line along the axis of length Zo, which is the ground state of 
unperturbed Hamiltonian. We now construct all possible higher excited states from this state following 
the tree structure outlined below. 

1. Start with an initial state (or, a path) with X coordinates varying from 0,1..., Zo and Y coordi¬ 
nate always zero. Length of the path is Iq. Assign a level index si and a serial number (say pJ 1 ) 
and. This is now a parent state at level s i. 

2 . Start from one end of the path p J 1 and consider each link of the path at a time. For each 
link, apply the appropriate BFACF move. Check if the move makes the path overlapping or 
backtracking. If not, the path created is a new path. Assign a particular serial number to each 
of them. 

3. repeat step 2 until one reaches, the other end of the path. 

4. In this step, a set of paths are created at a level, say 52 with length Zo + 2 which are say, q S2 
in number. All of these paths has the common parent (p^ 1 ). Assign the serial numbers to 
p{ 2 ,p 2 2 4 .. -P S 6 q S2 to each of these paths. 

5. Now, each of these paths p{ 2 ,p 2 2 • • -Pq 2 S2 will act as parents to create a set of paths in level 53 . 

6 . Start with path p^ 2 , repeat the steps 1 , 2 , 3 to create all possible paths of length Zo + 4. These 
are the paths at 53 with parent p \ 2 . 


10 
























7. Then consider path p^ 2 onwards and follow the step 1 and 2 to create a new path at level S 3 . For 
each new path created, check with already created paths at the same level, i.e S 3 , whether this 
path is indeed a new path or not. If yes, assign a new path number to it and if not, assign the 
path, with which it is overlapping, a new parent from which this new path has been created. 

8 . When exhausted with all paths at $ 2 , consider all the paths created in level S 3 as parent paths 
to create level S 4 following exactly same prescription as in 5, 6 , 7. 

9. Continue this tree structure to create states depending upon upto which order you want to 
continue in strong coupling correction 

The above algorithm is illustrated pictorially in figure [4] for a specific example, i.e for Z 0 = 3 and upto 
level S 3 



Figure 11 : Started with a string of length 3 units, all possible strings have been created upto length 
7 units. The strings of length 5 units are all created from the initial state. But for the next level of 
states, which are created from different parents of length 5 units can have multiple parents (denoted 
by more than one arrow meeting at a particular state). This feature remains valid for higher excited 
states which have not been illustrated explicitly. Note that, the arrows in the states are not relevant 
for SU( 2 ) string states, but are important SU(3) and higher rank groups 

The most important feature of the simulation algorithm is as follows: 

• First of all, the above algorithm is an exact simulation algorithm and hence is free from any 
statistical error. 
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• The simulation is in loop space and hence is free from any gauge redundancy. 

• starting from the strong coupling vacuum (or ground state), all possible loop configurations are 
built at each excited state contributing to correction of strong coupling result. Note that, this 
scheme involves no further truncation scheme unlike the other loop simulation techniques such 
as coupled cluster algorithm. 

• The algorithm creates states on an truly infinite lattice. 

The next part of the calculation involves evaluation of the perturbation expansion taking into 
account all the states simulated. Actually, for the sake of the present problem it is not necessary to 
store all the detail of the states as the non-zero matrix element of the interaction Hamiltonian is always 
the same as given in [22| The matrix element is non-zero only for the two states which are parent-child 
of each other. In the algorithm described above, we basically construct a parent-children directory for 
states and keep track of the excitation level of each. Finally, we utilize that to compute the expression 
given in © for perturbation correction order by order. The observations are discussed in the next 
section. 


5 Results 

First of all, starting with our effective Hamiltonian given in figure [4j we create all possible string 
states, starting from an initial string state of length Zo which is the ground state of the unperturbed 
Hamiltonian. We tabulate the number of string states produced at each energy level upto next 8 th 
excited levels in table [l] We plot this density of states in figure [5] and find out that the number of 

Table 1: Degenaracy factor at each excited level of unperturbed Hamiltonian 


length of the ground state (Zo) 

3 

4 

5 

No. of states with length Zo 

1 

1 

1 

No. of states with length Zo + 2 

6 

8 

10 

No. of states with length Z 0 + 4 

18 

30 

46 

No. of states with length Zo + 6 

64 

108 

176 

No. of states with length Zo + 8 

242 

416 

694 

No. of states with length Zo + 10 

896 

1584 

2720 

No. of states with length Zo + 12 

3290 

5938 

10448 

No. of states with length Zo + 14 

12134 

22192 

39732 

No. of states with length Zq + 16 

44884 

82932 

150474 


states do grow exponentially, the exponent giving notion of a non-zero Hagedorn temperature present 
in the theory m given by the values of the fitting parameter c in the last row of table [2] for different 
string states with different initial length (Zo = 3,4,5). 

Having generated all the string states connecting a pair of static quark-antiquarks, we calculate 
the perturbation expansion for the ground state energy of the system following (10) and © cor- 
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density of states vs energy level 



Figure 12: We fit the data in table [l] in function N = a + b e cn , where, N is the number of states at 
ttM excited level and a, 6, c are fitting parameters the values of which are given in table [ 2 J 

Table 2: Fitting N = a + b e cn , c denoting the Hagedorn temperature 


^0 

3 

4 

5 

a 

2.76577 

-4.04703 

-23.98411 

b 

1.28253 

2.18281 

3.55591 

c 

1.30787 

1.31815 

1.33164 


rented upto 16 th order taking care of contributions from all the states simulated (counting total of 
61535,113209,204301 for strings with minimum length 3,4,5 lattice units respectively). In figure |5j 
we plot the corrected ground state energy of the string state against length of the string state for 
various values of the coupling ?/, ranging from 1.5 to 1000. For each case we find a linear graph with 
the slopes denoting constant string tension. We further plot the slope of each straight line in figure 
[5] against the coupling for each case in figure [5| Extrapolating the graph to smaller values of the 
coupling, i.e, beyond the roughening transition point (around y = 1), we find the following estimate of 
the string tension as given in table [3] 

Table 3: Continuum extrapolation of string tension 


y 

T 

0.001 

0.688549 

0.0001 

0.670781 

0.000001 

0.645426 
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Strong coupling perturbation calculation of energy of a string of length L 



Figure 13: Slope of each graph gives the value of string tension for each case. 

6 Summary and Discussions 

This work, being the very first attempt exploit the prepotential formulation of lattice gauge theory 
from numerical perspective, contains an exact simulation algorithm for simulating all possible loop 
configurations of the theory and some preliminary results illustrating the viability of this new scheme 
of computations. We use BFACF algorithm to generate the loop configurations of the theory, starting 
from strong coupling ground state. Note that, each and every possible loop configurations are generated 
in this algorithm on an infinite lattice. Hence, except the strong coupling limit, no abrupt truncation 
of loop states is considered. This exact algorithm generates the loops staring from a parent state 
following a tree structure. The information of the tree structure or the network with all the loops 
acting as a node enables one to compute strong coupling perturbation expansion using the analytically 
calculated local dynamics of the theory in prepotential formulation. 

One of the most important feature of prepotential formulation is that, it gives local description of the 
theory as well as its dynamics. This enables one to split the link operator and hence the full plaquette 
operator in the Hamiltonian into a set of plaquette operators as given in figure |2j This splitting 
enables one to construct different effective Hamiltonians of the theory keeping all the symmetries of 
the system intact. Of course this choice of an effective Hamiltonian depends upon the problem one 
is interested in. For example, in this work we have chosen an effective Hamiltonian, which causes the 
dynamics of the theory only within the connected loops or strings. The disconnected loops do never 
appear in this modified effective theory. Moreover, our chosen set of plaquette operator acting on 
any loop state, always created a loop state with length two unit more or less. But there exists some 
terms in the Hamiltonian for which one plaquette action deforms a loop in such a way that the length 
remains same. We have particularly removed those terms from our effective Hamiltonian. Inclusion 
of such terms would have made the tree to grow sideways as well as opposed to the downward growth 
as illustrated in figure [4j Inclusion of such terms in the algorithm requires modification by including 
length changing BFACF moves together with appropriate change in calculational scheme. It is indeed 
important to move forward with this next level of complications by including those terms as well, in 
order to obtain the dynamics within all connected loops and to compare it with the results obtained in 
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Axial string tension vs coupling constant 



Figure 14: We fit the graph to the equation T = ao +a\ In y, with ao = 0.721492 and a± = 0.00550585. 

coupled cluster approximation. Moreover different effective Hamiltonians obtained within prepotential 
framework is worth investigating in order to establish the desired universality between these. Work is 
these directions are already in progress and will be reported soon. 

Finally, this algorithm is applicable beyond lattice gauge theory wherever loops are constructed on 
a square lattice. 
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